A GENERAL METHOD FOR DYNAMIC ANALYSIS OF STRUCTURES 
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A__Goneral Method for Dynamic Analysis of Structures 

I. Definition of the Ter m “Element* 

An element is defined as any structural unit whose degrees of freedom (DO PS) can be categorized 
as either interfac e DOFS or non-in terf ace (internal) DOFS. 

Interface DOFS are those DOFS through which the element is connected to one or more neigh- 
boring elements. Non-interface DOFS are internal to the element and do not directly couple to 
neighboring elements. 

The term “element” then, has a rather broad meaning. An element could be a fundamental 
structural unit such as a rod, a beam, a plate, etc. or it could be an entire structural component. 
Furthermore, the parameters of the element could be distributed or lumped. Figure l schematically 
illustrates the element concept. 

Any structural system can now be thought of as a composite of n such elements. The choice of 
elements is totally arbitrary and is a matter of user convenience. In particular, the user does not 
have to worry about the “size” of the element as is the case for example when using a standard 
finite element approach. This means that in general n is relatively small and little bookkeeping is 
necessary. 

Spatial periodicity of structures can be taken into account in a natural manner. These important 
advantages will be further clarified in the next few sections. 

II. Modeling of an Element 

Each element will be modeled using a set of assumed modes. In particular, a combination of 
interface constraint modes (ICM) and a set of interface restrained norm a l modes (IRNM) can be 
employed. Note that other types of assumed modes can be selected and should be investigated. The 
above choice is motivated by the Craig/Bampton approach to component modes synthesis and has 
several important advantages. 

I • Interface Co nstraint Modes (ICM) 

Let us assume the element has<j interface DOFS he ICM corresponding to DOF X^j 

is defined as the stati c deformation pattern of the element for X^y -/ and X Xi - 0 (for all L*j ). 
Note that many “shape functions” used in the finite element method are actually ICM. 

2. Interface R e strained Norm al Modes (IRNM) 

IRNM are the regular mass normalized mode shapes for the element with fixed interface DOFS. 

3 . Pis pi ace men t of a n Ele mer 1 1 

The displacement vector of a discrete element can be partitioned as: 



The non-interface displacement can be written as 
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wlu>re the first t.orm S represents a static deformation due to the interface displacements. The 
term X is best described as the difference 

~N ~ ~r 

Thc displacement vector X is now written as a linear combination of the IRNM 

x = & 4 

/V 


CM 


(») 


120 





whore N is a set of modal coordinates. 

The displacement vector X tan now be written as 

f-? 


ORIGINAL PAGE IS 
Of POOR QUALITY 


with 


X = 

[t 1 

1 4 

‘ X 


. s ni . 


LM 




( 5 ) 


(«) 


The columns of matrix "t clearly represent the ICM of the element. The actual general form of the 
matrix S^ x will be given shortly. 

The elements of the vectors and will be the generalized coordinates used to model a 
structural element. Equation (5) is al so valid fo r distributed parameter elements. 

Equation (5) indicates why the elements can be chosen with a large degree of freedom. If the 
matrices *6“ andi“*^ can be determined with enough accuracy then a legitimate element is found 
regardless of its size. For example entire beams, rods and plates can be considered as one element. 
Even large components could be considered as single elements as long as Eq. (5) can be adequately 
written. As a consequence, a considerable amount of bookkeeping can be avoided, thereby reducing 
the cost of the analysis. 

In the next section we will introduce the element equations of motion. 

II I . The Element Equatio ns pi Nipt jo n 

From here on we will work with discrete systems because relationships can be shown more 
explicitly. It should be kept in mind however, that all results are equally valid for distributed- 
parameter systems. 

The most general form of the element displacement vector X is: 


with 
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X can be written as: 
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Decause distinct interfaces by definition do not connect directly to each other, the element 
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The mass matrix has a similar form. A notliQF f<JPflQBrr6yAWiTiX to vector x in Eq. (9) is: 
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with partitioning corresponding to Eq. (9) 
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The form of Eq. (10) suggests that 
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where are IRNM and the partitions in Eq. (14) correspond to Eq. (13). 

Using Eqs. (II), (12) and (14) the mass and stiffness matrices for an element corresponding to 


coordinates x 7 are: 
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IRNM ortliononnality conditions 
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The approach followed to model the element reveals a unique perspective on a serious shortcom- 
ing of the standard finite element techniques when used to construct element mass matrices. Indeed, 
it is clear that a standard finite element mass matrix (consistent) represents a Guyan reduction in 
which the internal degrees of freedom (non-interface) are eliminated, i.e. all coordinates are 

neglected. As will be seen later, such an approximation is valid only when the elements art* '‘small*' 
enough so that indeed the coordinates ^ have no effect on the response of the structure. 

In the next section we will discuss the formation of the system equations. 
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V V . Tho_Systern Equations 
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In this section we will demonstrate how elements arc assembled into global components. As an 
example, consider the three elements in Figure 2* Equations (16-20) represent the general form of 
the element equations. The displacement vector for clement L can be written as 


Kl j 

CJ. 
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In general, the vector represents a set of loca l displacement coordinates. In order to enforce 
displacement compatibility between the element interfaces it is necessary to transform into a set 
of global coordinates ^ x t ^ ^ ^ 
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where is a geometric transformation matrix. Note that no transformation is necessary for 
because elements only connect through interfaces. Using equations of the type Eq. (24) \ 


rewrite Eqs. (16-20) in the following form 
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whereM^M^and are easily found. 

The first step in the assembly process is to connect clement l to element 2 by requiring that 
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X * 


for all times 


(27) 


The unc oupled equations of motion for the 1-2 component can be written as follows 
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where the damping and forcing termsTare omitted since we are only interested in the eigenvalue 
problem. Note that for element 2 we have two interfaces and we denoted 
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Taking into account Eq. (27) wc ran write tin* coupled equations of motion for component 1-2 as 
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At this point, element 3 should be connected to the 1-2 component. Before doing this. Equation 
(30) will be manipulated into a form identical to Equations (16-17). This is an important feature of 
the presented technique. Indeed, it will allow us to consider the 1-2 component as a standard element 
as defined in Section I. In other words, the resulting equations will represent the dynamics of element 
1-2. In addition, the equations will be in a form suitable to apply an escalator type eigensolver to be 
discussed later. 

2. Transformation of the Interface Displaceme n t Vector 

The first step in reducing the 1-2 component into a single element is to introduce the following 
transformation 
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This transformation is entirely equivalent to Eq. (2). The term S, TX 3^ rep resents that part 
of which is due to the displacement of the interface 12. The term 5?^ ' s displacement of J 
with respect to the interface 12 as seen by an absolute observer. Note that the matrix^VA^is never 


singular because it represents the stifFness matrix of a rigid body restrained system. 

The displacement vector J £ can be written as a linear combination of IRNM as follows 
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where the eigenvalue problem solved is 
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can now be substituted into Kq. (30). yieldin'; 


tu lt 

*3/ 




H 

X 


y». 


32 




42 


**>3 

I 

o 


tv\ 


14 


M. 


14 


0 

I 


£ 

t, 

ut j 




o 

o 

0 


I 

0 

* 1 , 

0 

0 


] 

< 







J 


(33) 


(3.1) 


(33) 


0 

0 

0 


0 

0 

0 


rrZ 


c 

I'r 


where ... ran ,M> oas,| y found. 


% 

i£ j 


(:«i) 


125 



3. Thc JIlNM Kigcnvaliic Problem 

The second step in reducing Kq. (:«») to a form identical to K«|s. (Ki-17) is the solution of the 
eigenvalue problem 
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Using the transformation 



(37) 


(38) 


(39) 


with 


, P,2~ *>/«• »n 

K - 


(»«) 


Equation (39) is now equivalent in form to Eqs. (16-17). It essentially represents the equations of 
motion for the 1-2 component viewed as one single element. If along the way no modes are truncated 
then Eq. (39) involves no unusual approximations. In particular, for the continuous equivalent, flu* 
equations corresponding to Eq. (39) are still “exact”. 

The eigenvalue problem (37) has a very special form and an optimum solution will bo discussed 
in Section V. 

An important aspect of the present technique is the truncation of the mode sot according to a 
preset frequency. This preset frequency must reflect the frequency content in an element necessary lo 
obtain the desired fidelity in the overall system model. How the “element” cut-off frequency compares 
to the “system” cut-ofT frequency is still a matter of research or “experience". 

Once the appropriate truncation is performed, a third clement can he added through interface 
(2. Note that the number of degrees of freedom is already reduced. The eigenvalue problem (37) is 
small and can be solved very efficiently as will be shown. 

It is also important to note that because of the transformation (31) we do not loose any accuracy 
in the rigid body and static properties of the system when modes are truncated. 

I. Coupling of Element 1-2 and Moment .3. 

T3ie next step is to couple element .3 to element 1-2. This process is very similar to the one 
already described. From Eqs. (25-26) we can obtain for C 9 S 



The uncoupled equations of motion for the 12-3 component can be written as 
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where we used Eqs. (39) and (41). Note that the interface 12 of element 3 is part of the model if no 
other elements are to be connected. However, if other elements are to be connected, this interface 
must be made explicit. 

We now must impose the compatibility condition 



Because this is the last element, no further transformations of the interface displacement vector is 
necessary. Using Eq. (34) we can write 

•£„+«?, M*. M,\l fgM K„+K,f o o 

X 0 i 7f U 0 CS X 0 

Ml ° 1 JlSj [ « 0 <S\ 

which represent the coupled equations for elements 1-2 and 3. 

Next, we can solve the small eigenvalue problem. 
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so that „ ^ 

* ± 2 X ( «o) 

Note that no “bar” is necessary over the above quantities, because if the system is free, the rigid 
body properties will be incorporated in Eqs. (15) and (46). 




Finally, note the simple form of this system of equation (47). First, generally truncation will redui e 
the size significantly. In addition, the specific form will allow lor application of the escalator algorii lirii 
as discussed in Section V. This eigenvalue solver will yield a set of system frequencies jfcj^tud svstem 
modes ^5" with a minimum of effort and cost. 

Ilefore we discuss the special eigenvalue solver, remarks are in order: 

(I) So far, we described how the elements are coupled together. It is a matter of repealing tl ie 
same procedure for each added element. Fiach time truncation is used on the element level as 
well an on the level of the current system. A series of relatively small eigenvalue problems is 
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solved using a very efficient eigenvalue problem solver. The entire system is gradually built 
up, keeping only the necessary frequency fidelity at each stage. Several interesting questions 
arise, for example: (a) what is the ideal tradc-ofF between element frequency truncation and 
current system frequency truncation, and how do both relate to the overall desired fidelity? (b) 
Which is the optimum way of solving the successive eigenvalue problems? Should we wait until 
several elements are collected before a current system eigenvalue problem is solved? This is an 
important question since it affects storage, cost and accuracy. 

(2) The manner in which the elements are coupled together makes this technique ideally suited to 
handle spatial periodicity. Indeed, Eq. (39) shows that each current system can be considered as 
an “element”. If care is taken, it is possible to use the same element over again, without adding 
significantly to the cost. For example, a periodic truss can be started with one clement which 
is truncated according to a desired frequency. This model can now be doubled and truncated 
again. This truncated 2-element model can now be doubled again, to yield a -(-element truss, 
etc. 


V. The Eigenvalue Problem Solver 


In this section we will describe an eigenvalue problem solution technique which is particularly 
suited for our purposes. 

First, let us consider the following special eigenvalue problem: 
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(19) 


This problem has a diagonal stiffness matrix and a unit mass matrix except for the first, row and 
column. The Ml values are such that the mass matrix is positive definite. 

The characteristic equation of this problem can be written as follows 
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(50) 





In other words, we have isolated the eigenvalues of the system represented by Eq. (19). Note that 
property (st) again shows that for^-aik* ^.becomes a root of the system. Property {ft) allows us to use 
for example the Newton- If aphson toennique to find the actual eigenvalues^ ' . This iteration scheme 
converges quadratically provided a good initial value is found. Without going into detail, at this 
point we can say that property If/) allows for a very accurate initial value for each of I lie eigenvalues 
. Therefore, convergence is extremely fast, often after two or three iterations. 

Once the eigenvalues are determined, it is computationally a trivial mailer to obtain the cor- 
responding eigenvectors. Also, multiple roots are no problem. This algorithm was programmed but 
not yet optimized. Nevertheless, on a VAX/ 780, it requires only 10 seconds of CIM' time to solve 
a problem with n 150, which is a size far beyond our needs. As part of the presented dynamic 
analysis technique this eigenvalue problem must be solved for each interface dot. 


VI. Su mrnary and Cone Ijisions 

The presented research deals with the development of a dynamic analysis method for structural 
systems. The modeling approach is essentially a finite element method in the sense that the structure 
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is divided into n “elements”. An “element” is defined as any structural unit whose dofs. can be 
categorized as either interface or non-interface dofs. The term “element’' then, has a rather broad 
meaning. An element could be a fundamental unit such as a rod, a beam, a plate etc. or it could be 
an entire structural component. Furthermore, the parameters for the element could be distributed or 
lumped. The choice of elements is totally arbitrary and is a matter of user convenience. In particular, 
issues of accuracy and convergence do not enter on the level of element choice as is the case in a 
standard finite element method. This means for example that bookkeeping is reduced to a minimum. 

Each element will be modeled using a set of interface constraint modes (ICM) combined with a 
set of interface restrained normal modes (IRNM). If the ICM and IRNM can be found with enough 
accuracy, then a legitimate “element” is defined. For example, entire beams, rods and plates can be 
considered as one element. Entire components can be made into one element in an off-line manner. 
Moreover, the entire system can be modeled as if it was represented by partial differential equations. 
Note that standard finite element techniques in general only use ICM ( - shape functions) which 
leads to very interesting and useful insights into important shortcomings of these techniques. In 
particular, the problem of an accurate mass distribution is addressed by this new approach. 

The element assembly process is essentially the same as in the standard finite element method. 
However, each combination of elements is automatically converted into a single element. This pro- 
cedure is based on static condensation without loss of accuracy. This feature is very important and 
allows for each structural unit to be interpreted as an “element”. It also allows for the stiffness 
matrix to remain diagonal. 

The next step is the solution of the system eigenvalue problem. The procedure calls for the 
sequential solution of a number of small eigenvalue problems based on a truncation principle for 
IRNM. In addition, the form of these eigenvalue problems is very simple such that an escalator type 
of eigenvalue problem solver can be used which is extremely cost-effective and fast. The response, 
loads, etc. calculations are rather standard, but also benefit from the approach in terms of accuracy 
and cost-effectiveness. The groundwork for this technique is in place and is currently supported by 
the AFWL/ARBH Kirtland AFB, NM. Some of the advantages of the new technique are: (l) The 
problem of Order Reduction is believed to be solved. The technique implicitly reduces the system 
order. Whenever an element is added only information necessary to obtain a prescribed fidelity in the 
system model is retained. (2) Very accurate. In fact, if desired, “exact” solutions in the distributed 
parameter sense can be obtained for any structure. (3) Fast and cost-effective. This is due to the 
small number of elements: the solution of a series of small eigenvalue problems instead of one large 
problem; the special nature of these small eigenvalue problems combined with the cost-effectiveness 
of the escalator eigenvalue problem solver. ( I) Applicable in general. In particular, extremely large 
structures do not pose a problem. Once a model is agreed upon, any number of modes and frequencies 
with any degree of accuracy can be computed. (5) It is anticipated that Micro Computers can be used 
to solve even the largest of problems. This is due to the small bookkeeping effort and the sequential 
nature of the solution. (6) Spatial periodicity ran be taken into account in a natural manner. (<) 
It is anticipated that this method will be useful in areas like control optimization, identification 
and possibly non-linear phenomena. The feasibility of this technique as well as several of the above 
advantages have been demonstrated with several examples. 
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